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Abstract 

Within the finite-temperature Greens-function formahsm we study the equation of 
state of a hot interacting pion gas at zero chemical potential. Employing real- 
istic vrvr meson-exchange interactions we selfconsistently calculate the in-medium 
single-pion selfenergy and the vrvr scattering amplitude in the quasiparticle approx- 
imation. These quantities are then used to evaluate the thermodynamic potential, 
il7r(T), from which the state variables: pressure-, entropy- and energy-density can 
be derived. In contrast to earlier calculations based on the low-energy Weinberg 
Lagrangian we find an overall increase as compared to the free-gas results. We also 
consider the possibility of a dropping p- meson mass as suggested by the 'Brown- Rho 
Scaling' law. 
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1 Introduction 



In recent years much effort has been put into the study of ultrarelativistic heavy-ion 
collisions (URHIC's). The main objective of these experiments is to create a new state 
of hadronic matter, the Quark-Gluon plasma (QGP). However, it is not settled if the 
current generation of experiments at the BNL-AGS and CERN-SpS is able to generate 
sufficiently high energy densities to make the transition. On the other hand, lattice gauge 
calculations, although not yet at the stage of making accurate predictions, suggest the 
occurrence of a different kind of phase transition associated with the restoration of chiral 
symmetry. This 'chiral phase transition' is expected to occur prior to the deconfinement 
transition at a critical temperature of = 140 — 170 MeV. Precursors of chiral symmetry 
restoration may well establish before that, e.g. the dropping of vector meson masses as 
proposed by Brown and Rho [|l], 0] based on chiral effective Lagrangians and further 
corroborated by Hatsuda et al. within the QCD sum rule approach. 

In central collisions at the AGS, the midrapidity region is characterized by highly 
compressed nuclear matter due to an almost entire stopping of the Lorentz-contracted 
colliding nuclei. Even at SpS energies, where the stopping power is much less, one en- 
counters sizable baryon densities at central rapidity, leading to an appreciable impact on 
the in-medium properties of the produced secondaries (mostly pions). At the Brookhaven 
Relativistic Heavy Ion Collider (RHIC), on the other hand, baryon densities in the cen- 
tral zone are expected to be very small, such that, after hadronization, this zone will be 
populated by a dense pion gas. 

Many aspects of a hot, interacting pion gas have been studied in the literature so far, 
e.g. single-pion 'optical' potentials 0, ^, ||, 0], mean free paths f^, in-medium rni cross 
sections |10|, H , the equation of state [Tl|, |T2| , hydrodynamic properties [l^ as well 



as numerical solutions of a bosonic Boltzmann equation [|T3] simulating the dynamics of 
URHIC's prior to freeze out. 

In this article we want to concentrate on the equilibrium properties of a thermal pion 
gas, i.e. the equation of state (EOS). Within the relativistic virial expansion, restricted to 



two-body collisions, Welke et al. [11| have previously employed an empirical parametriza- 
tion of the experimental vacuum vrvr scattering phase shifts in s- and p-wave as their 
basic input for calculating the number-, energy-, pressure- and entropy-densities of an 
interacting gas of pions. They have found that the net effect of the tttt interaction stems 
from the resonant JI=ll-(p-) channel leading to an increase of the thermodynamic state 
variables, relative to the free gas, for temperatures T~100 MeV. Quantitatively this con- 
tribution is comparable to what one would expect from an admixture of free p mesons. In 
other words, the rather sharp p resonance in the tttt interaction resembles the contribution 



from free p mesons, a mechanism well-known from the Beth-Uhlenbeck formalism |]16 



Bunatian and Kampfer [0] have pursued a rather different approach, based on finite 
temperature Green's functions in the imaginary time formalism. As the tttt interaction 



they employed the Weinberg Lagrangian ||T^, known to account for low-energy s-wave vrvr 
scattering. As a result of their selfconsistent calculations in the Hartree approximation 
the thermodynamic state variables exhibit a decrease for temperatures T~ 150 MeV. 
The aim of the present paper is to recalculate the EOS by both employing a realistic tttt 
model capable of describing the vacuum scattering data up to rather high energies and 
by taking into account the in-medium modifications of single-pions as well as of the vrvr 
scattering amplitude selfconsistently and to all orders. Our formalism for calculating the 
thermodynamics will be the one used by Bunatian and Kampfer, but extended beyond 
the Hartree approximation. The vrvr interaction we use is the meson-exchange model de- 



veloped by the Jiilich group [|T8[ which, on a microscopic level, describes the tttt scattering 
data in the various partial waves up to CMS energies of 1.5 GeV. 

Our article is organized as follows: in sect. 2 we present the formalism for calculating 
the thermodynamic potential f2^(T). In sect. 3 we derive from it the thermodynamic 
state variables, i.e. the pressure-, entropy- and energy-density, and discuss the numerical 
results within the Jiilich vrvr model with no in-medium modifications applied to the meson- 
exchange potentials (Brueckner theory). In sect. 4 we consider the possibility of additional 
medium effects on the exchanged mesons in terms of a dropping p-meson mass. In sect. 5 
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we summarize and make some concluding remarks. 



2 The Thermodynamic Potential of Interacting Pi- 



ons 



The starting point of our analysis is the well-known expression for the thermodynamic 



potential of a gas of interacting bosons at finite temperature |]T9[. In case of interacting 
pions it reads 

n^iT) = n^{T) + n^^iT) (i) 

where the first term represents the quasiparticle contribution while the second term arises 
from interactions between the quasiparticles. For thermodynamic consistency both terms 
have to be considered. 
In terms of the single-pion propagator 

D^{u+,k) = (^ul - ml - k'^ -J:^{uj+,k)j , (2) 

(a;+ = u + irj), where T,.„{uj+, k) denotes the pion selfenergy, the quasiparticle contribution 
is given by 

= / ^ / ^ r(^) Im{ln[-D;iK, k)] - D.K, fc)S.K, k)] (3) 

where f""{uj) = {exp{uj/T) — 1)^^ is the thermal Bose factor. The interaction contribution, 
Q^^[T), is the sum of all 'skeleton diagrams' arising from a perturbative expansion of the 
scattering amplitude. Within the Matsubara formalism the contribution from n-th order 
meson exchange can be written as 

Q 1 r rl^k 

^^^HT) = ^^i-T) J ^,i:DA^u,k) . (4) 

This can be pictured as closing the two external legs of the pion selfenergy. The factor 
l/4n corrects for overcounting |T^ induced by the possible ways of regenerating the n-th 



order pion selfenergy by cutting a single-pion line in the diagrams of r2^^(T). For the case 



of TTTT contact interactions as where considered in refs. [0, ^ the appropriate counting 

factor is l/2n due to the absence of u-channel exchange graphs. 

Let us first turn to the calculation of the two-body interaction contribution f2^7r(T). 



2.1 Lowest-Order Contributions to ^}-J^^^{T) 

To lowest order [n = 1) Qt^t^{T) has already been calculated in ref. [|1^]. In the context of 
a vrvr meson-exchange potential it reads 

f ^r(co')ImD^ico',p)Mi'Jico + co',k,p) , (5) 

J 71 

where M^]^ denotes the spin-isospin averaged, first-order forward scattering amplitude. 
Employing the quasiparticle approximation (QPA) for the single-pion propagator, 

TT 

ImD^{uj, k) = - — [S{uj - Ck) - 5{uj + Ck)] , (6) 
eq. (5) can be simplified as 



oj 27T^ J (27r)2 2ek 2ep 
{MW(e, + Ck, k,p)[r{ek)r{e,) + {r{ek) + l){r{e,) + 1)] 
+Mi'J{ek - e„ k,p)[nek){ne,) + 1) + (r(e,) + l)r(e,)]} , (7) 

where x = cos 6 ,6 = l{k,p). The quasiparticle energies are determined from the 
selfconsistent solution of the Dyson equation: 

el = ml + k'^ + ReJ:^{ek,k) . (8) 

The right-hand side (rhs) of eq. (7) diverges due to the various 'I's appearing in the 
occupation factors. However, when the thermodynamic potential is defined relative to 
the physical vacuum by only considering the difference between the values at T = and 
finite T then it remains finite. This is achieved by removing the 'I's in eq. (7) [0. Thus 



3 r edk r p^dp 
^ si 2^7 M^^^ ^"'^^^ 

-1 

1 rfx{M«(ep + efc, + M«(efc - e^, . (9) 

+1 

In the following we shall suppress the tilde and always refer to subtracted quantities. 
For practical evaluations, the M-amplitude has to be transformed into the CMS of the 
two scattered pions, thereby neglecting the medium dependence of M^t^ on the total 
momentum P = k + p oi the pion pair. Then the total CMS energy is identified as 

E^ = s = {ek±epf-{k + p)^ (10) 

and, by requiring Lorentz invariance, the relative CMS momentum of the pions is obtained 
as 

= l[\{s - {el - e) - {el - p')f - {el - fc^)(ej - p')] . (11) 
Using ds = —2pkdx we can rewrite eq. (9) as 

with 

^max 

I^A/B = / EdEMi'J{E,q,q), 

rrnn 

ET = {e,±e,r-e-p' + 2kp, 

E^l,^ = {ek±epf-k^-p^-2kp, (13) 
which is the main result of this subsection. 



2.2 Higher-Order Contributions to r27r7r(T) 

Within the meson-exchange framework, the full Mn-Tr^amplitude is obtained by solving 
a Lippmann-Schwinger type equation (LSE), which, in our case [^, results from the 
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Blankenbecler-Sugar reduction |21| of the covariant Bethe-Salpeter equation. In a given 



partial- wave/isosp in channel J I the scattering equation reads 

Mii{E,q,,q,) = V^'^ {E,q,,q,) + Vi'{E,q„q) G^E^q) Mii{E,q,q,) ,(14) 

^ ' 

where (7^,^ denotes the two-pion propagator of the intermediate state and Kttt = 

are the Born amplitudes (pseudopotentials) derived from an effective meson La- 
grangian. Schematically, the scattering equation (14) can also be written as a perturbative 
expansion 

M = V+VGM 

= V +VGV + VGVGV + ... 

= M^^) + M(2) + M^^) + . . . . (15) 

Beyond the Born approximation, i.e. for n > 2, the M-amplitude acquires an imaginary 
part due to the intermediate two-pion states being on-shell. Consequently, S^"^ also be- 
comes a complex quantity. In order to evaluate eq. (4) for n > 2, we therefore employ the 
standard procedure of inserting the spectral representations of the single-pion propagator 
and the selfenergy: 



duj' ImDT,{iu', k) 
n Zy — uj' + if] 



^^:\z.,k) = -/^ '"^^-"^("'.'^ -i?eE(")(0,fc) 

J 7T Z^ — UJ + 17] 



( -Im^^:\u, k){ + ^) . (16) 

J 7T Z,, — UO + in UO 



According to the findings in ref. |^ the dispersion relation for the pion selfenergy has to 
be supplemented with a subtraction at zero energy. 

The Matsubara sum in eq. (4) can now be performed analytically. The imaginary part of 
the pion selfenergy has been derived in ref. p2[ and we can immediately generalize the 
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result to the n-th order contribution: 



/mEW(^+,^) = -| -^ImMil\u;+W,k,p)ImD,{u;'^,p)[r{u;')-r{u;+u;')] . 

(17) 

Injecting eqs. (16), (17) into eq. (4) results in the following, exact expression for the n-th 
order scattering contribution to the thermodynamic potential: 

^(r,\ 3 1/" d^k f d^p f dE' f du f du' ^ ^ , . ^ ^ ^ , , , 

ImM^^ iE,k,p){ iu;^+u;'-E')iE'-u;') ^ " ^^^^ 

Note that the imaginary part of the rhs , given by the 5-function part of 1/ (a;+ + a;' — E') , 
vanishes due to the simultaneous disappearance of the second occupation factor - as it 
should, since flT^^{T) is a real quantity. Applying the QPA, eq. (6), and performing the 
trivial angular integrations leads to 

2[F^{E',k,p) + G-{E',k,p)] , (19) 



where 



F''{E',k,p) = 



G-{E'Xp) = 



[r{e,) - riE')][e,rie,) - {E' - e,)r{E' - e,)] 
{ek + ep-E'){E'-ep) 
^ [r{e,) + r{E')][ekr{ek) - {E' + e,)r{E' + e,)] 
{ek + ep + E'){E' + ep) 
[rjep) - r{E%e,r(e,) - {E' - ep)r{E' - e^)] 

{-ek + ep- E'){E' -ep) 
irjep) + nE')][e,neu) - {E' + ep)r{E' + e,)] 

(e,-ep-E'){E' + ep) " ^ ' 



The second term in F^ arises from the negative-energy part of the £^'-integration. It is 
obtained from the first term when replacing E' {—E') and introducing an overall minus 
sign due to the symmetry ImMT^T^{—E' -\-ir]) — —ImMT^T^{E'-\-ir]). Similarly, the function 
is generated from F'^ when replacing — > {—ek) with an overall minus sign due 



to the antisymmetry of the single-pion spectral function, eq. (6). Since the 4-momenta 
(cu, k), (uj'jP) enter eq. (18) on equal footing we obtain an additional factor of 2 for the 
negative energy contributions of ImD-,^{(jj\p). 

It is not to be expected that, evaluating Vt^^}{T) order by order, leads to a convergent 
series. For the LSE (15) the Born series does not converge and, although there appears 
an additional factor 1/n when summing the thermodjTiamic potential 

oo 

f]..(r) = Efi(';)(T), (21) 

n=l 

it may not change the convergence behavior of the series. In particular, the s-channel 
pole graphs (i.e. genuine vrvr-resonances) turn out to be problematic, since they exhibit 
non-integrable singularities in each order n > 2. In the following we will discuss them in 
more detail. 

2.2.1 s— Channel Pole graphs 

In the Jiilich model of Lohse et al. [0] a genuine resonance is characterized by a separable 
Born term of the form 

V^'iE', gi, q^) = v^^MDl{E')v.U(l2) (22) 
with the bare resonance propagator 

Dl{E') = ^—r^ , (23) 

m(j°) being the bare mass of the resonance a with spin-isospin JI. The vertex functions 
contain coupling constants, isospin and form factors as well as kinematic factors specific 
to the spin-momentum part of the coupling. 

As is well known, the LSE can be solved analytically for the case of a separable poten- 
tial which yields 



oo J f rpl\ 

n=i E'^ - [mh'Y 



£;'2-(mL°^)2-4(E') 



(24) 



where 

UE') = J qHq v^Uq? G^^E', q) (25) 

is just the intermediate tttt bubble. In eq. (24) the summation index is chosen such 
that n — 1 corresponds to the lowest-order (Born) term. Apparently the n-th order 
contribution in eq. (24) exhibits a pole of order n at £" = m^^^ , leading to divergencies in 
the ^^'-integration of eq. (19). The infinite sum of tttt bubbles, on the other hand, simply 
renormalizes the bare propagator generating a mass shift and a finite width, given by 
ReIa{E') and ImIa{E'), respectively. As we shall show now, this implies the possibihty 
of summing the VL^^'}^{T) contributions to all orders, thereby avoiding any divergencies. 
The evaluation of eq. (21) for pure resonance scattering involves the following sum: 



n=2^ 



OO 



= (26) 

where we have added a zero in form of the {n — l)-term (Born term), which does not 
generate any imaginary part. Using the identity 



OO 



ln(l-x) = -ETr' (27) 



n=l ^ 



we are now able to perform the infinite sum in eq. (26). The resulting thermodynamic 
potential reads 



^ 3 /■ k'^dk f p^dp /, JdE'^. . . . 



2x 



■1 



with 



^"^^TIW^^^- J'^^^\)J 2 [F-{E',Kp) + G-{E\k,p)] , (28) 
-'al-C/ j E'^ — {rria ) 



Sna)=^.2i±2M±21 (29) 
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being the spin-isospin weighting factor of the corresponding tttt resonance channel. When 
numerically evaluating eq. (28) in the Jiilich Model, one obtains results far too large to 
be realistic. A closer inspection shows that, for fixed momentum k, the p-integration 
dominantly picks up high-momentum components. The reason for that can be traced 
back to the weak suppression in the (dipole) form factors entering the vertex functions 

^7r7ra(?) ; 

In the limit of large p-values, p':$> k, the CMS momentum q behaves like oc p, i.e. one 
essentially looses four powers of momentum in VT^T^ai^Y oc F{q)^. Note that the occupation 
factors F'" and C", eq. (20), contain terms independent of /'^(cp) such that no 'thermal' 
suppression is at work. To overcome this unrealistic behavior we have decided to replace 
the form factor at each vertex by 

.(.)^;n^mri^(|f±£|)(|f±£|). (31) 



which restricts the relevant momentum ranges in eq. (28) to a realistic domain of p. A; ~ Aq, 
(e.g. Ap = 3.3 GeV in the Jiilich model). 

2.2.2 t— Channel Vector-Meson Exchange 

In the Jiilich model the dominant contribution to the low-energy tttt interaction is gener- 
ated by t-channel exchange of the p meson. The corresponding Born amplitude is of the 
non-separable form 

VfJ^{E\ gi, q^) = v,^,,,{E', q,, q^) D,{E\ q,, q^) , (32) 

where 

D,{E',q^,q2)^-^ (33) 
t — nip 

is the effective p t-channcl propagator with 'renormalized' (physical) mass, mp=770 MeV, 
and 4-momentum transfer t — {q2 — qi)^- The non-separability evades both an analytical 
solution of the scattering eq. (14) and a closed summation of the thermodynamic potential 
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eq. (21). Thus we are forced to evaluate the t-channel exchange contributions to O^^(T) 
order by order. Fortunately there are two circumstances that will allow us to get a realistic 
estimate: 

- the p t-channel exchange propagator, eq. (33), does not exhibit any poles in the 
kinematical region relevant for calculating the n-th order contribution to VL^^T^iT), 
thus ensuring the absence of divergencies order by order; 

- as it will turn out from the numerical results in sect. 3.3 the series eq. (21) rapidly 
converges when summing the closed ladder sum (e.g. the third order contribution is 
typically 10% of the second order contribution). 

When evaluating eq. (19) by using the expansion in eq. (15) with the pseudopotential (32) 
we encounter the same problem as was pointed out at the end of sect. 2.2.1, namely 
unrealistically large contributions from very high momenta. Therefore, as in the case of 
pole graphs, we replace the form factor entering the vertex functions in eq. (32) by 

F{q,qf^F{k)F{p) . (34) 

By a slight readjustment of the form factor we made sure that the exact n-th order p 
exchange amplitude is reasonably well reproduced when using this separable approxima- 
tion. 

Another characteristic feature of the Jiilich model is the occurrence of a scalar-isoscalar 
bound state in the coupled KK channel, just below the KK threshold. This bound state 
is generated by a strong attraction between kaon and antikaon, mediated by t-channel 
exchange of the vector mesons p(770), a;(782) and 0(1020). It leads to the characteris- 
tic sharp rise in the tttt —>■ in: phase shifts 5°° in vicinity of the KK threshold and is 
interpreted as the /o(980). To incorporate this resonance-hke feature of the tttt phase 
shifts in our calculation of Q,-^t^(T) we proceed as follows: rather than summing the ladder 
diagrams involving the KK intermediate states and their direct interaction we simulate 
the /o(980) as a genuine resonance, i.e. we choose the resonance parameters such that 
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the full 5°° phase shifts can be described without any coupling to the KK channel. The 
fully iterated /o(980)-polegraph contribution to VLt,t,{T) is then calculated as outlined in 
sect. 2.2.1 (eq. (28)). 

Before ending this section we should point out that we have neglected any interfer- 
ence terms between the s-channel pole graphs ('bubble' sum) and the t-channel exchange 
graphs ('ladder' sum). The following two arguments make us believe that this is a rea- 
sonable approximation: 

- the TTTT interaction in the resonant channels (i.e. the JJ=11 channel with the p(770)- 
resonance and the J/=20 channel with the /2(1270)-resonance) is largely driven 
by polegraph contributions such that the additional incorporation of t-channel ex- 
change processes is expected to have a quantitatively small effect; in the case of the 
p- and /2-channels the p t-channel exchange is attractive and leads to a rather small 
change of the resonance when evaluated with the combined ladder-/bubble-sum. 
Our results for the polegraph contribution, fi^^^(T), should therefore be considered 
as a lower estimate of the exact contribution; the situation is less clear for the (sim- 
ulated) /o(980) in the J/=00 channel, but as it will turn out from the numerical 
results the /o-bubble sum gives a very small contribution by itself; 

- as already mentioned in sect. 2.2.2 the p t-channel exchange contribution to VLt^t^{T) 
converges rapidly with the (r;,+l)-th order being typically down by 10% from the 
n-th order. 

We now turn to the evaluation of the single-(quasi-)pion part, f2^(T), of the thermody- 
namic potential. 



2.3 Quasiparticle Contributions: Q.^{T) 

In evaluating the contributions to Ql^^iT) we shall also apply the QPA to the medium- 
modified single-pion spectrum, given by eq. (2). We essentially follow the steps outlined 
in ref . p4 . 
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First we make use of the analytic properties of the one-pion propagator to transform the 
negative energy part of the (^-integration to positive values, namely 

D^{uj^,k) ^ Dl{uj-,k) . (35) 

Using the identity 

Im[\n{z)] = -/m[ln(z*)] (36) 

for any complex number z, and removing the infinite vacuum part (as discussed in 
sect. 2.1) results in 

f^^(r) = / 7^ / V 2^^^^^^ Im{H-D-\uj^, k)] - D^iu;^, fc)E.(^+, k)} . (37) 







The first term in the braces can be rewritten by using the identity 

Im[ln{-D;\u, k))] = 7re[ReD;\u, k)] - arctan[^^^p^] . (38) 

Neglecting the imaginary part of (which we will call scheme 'a') eq. (37) can be 
simplified as 

Q^'"(r) - ^ / dke{T\n{l - e-'^/^) - -L/-(e,)it:eE,(efe, k)} . (39) 

The first term corresponds to the contribution of 'free' quasiparticles with modified dis- 
persion relation — e{k), whereas the second term arises from an additional, dynamically 
generated, mean field. Their qualitative behavior is in line with the naive expectation 
that the pressure, p — —D,, is reduced for heavier particles but increases for a repulsive 
mean field and vice versa. 

One can go one step further by including the imaginary part of the pion selfcncrgy in 
eq. (37) (which we will call scheme 'b'). Then fl^{T) receives further contributions and 
reads 



27r2 



Ukk^ f^riu) ^ ; • (^o) 

J -l ^ ^ '{u;^-el)^ + {ImE,{ek,k)y ^ ' 
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We shall consider both approximations, eqs. (39) and (40), in our numerical calculations 
discussed in the next chapter. 



3 Thermodynamic State Variables and Numerical Re- 
sults 

3.1 Pressure, Entropy and Energy 

Having the thermodynamic potential at hand, the state variables: pressure-density, p{T), 
entropy-density, s(T), and energy-density, e(T), can be calculated by means of the stan- 
dard thermodynamic relations at zero chemical potential: 

p{T) = ~n{T) (41) 
siT) ^ -'^ (42) 

e(T) = TsiT)-p{T). (43) 

As elaborated in refs. [^, |19| the partial derivative w.r.t. T in eq. (42) only acts on 
the explicit temperature dependence of the occupation factors due to the stationarity 
condition 0^9 



Thus we obtain for the quasiparticle contributions to the entropy from eq. (39): 

s^'^T) = ^ / dke{{l + r(efc)) ln(l + r(e,)) - /"(e^) ln(r(e,)) 

-^^^^eS.(e.,A:)}, (45) 
or, when including a finite width for the pions, from eq. (40): 

s'^'\T) = / dk^iil + r(e,)) ln(l + r(e,)) - r(e,) ln(r(e,)) 



27r2 



arctan 



71 or uj'^ — ujI — ReJ2.„{ek, k)' 
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oo 







, duj df^{uj) ImJ:^{ek,k){uj^ -ujD 
' 71 dT (cu2-e2)2 + (/mS^(e,,A;))2^ ' ^ ' 



which is consistent with the expressions quoted in ref. [^]. In our framework of meson- 
exchange interactions, the skeleton contributions to the entropy (denoted by —d^/dT in 
ref. arise in lowest order from eq. (12), 



for the t-channel p exchange graphs from eq. (19), 

2[ + ] , (48) 

and an analogous expression for the s-channel pole graphs from eq. (28) for s,r7ra(7'). 
To check our results with respect to thermodynamic consistency, we will also consider an 
alternative way of evaluating the pressure density, namely 

T 

p^{T) = J dT' s^{T') . (49) 



3.2 Selfconsistent Brueckner Scheme 

As pointed out in sect. 2.1 the in-medium pion dispersion relation (or pion selfenergy) has 
to be determined selfconsistently from the Dyson eq. (8). In the QPA the pion selfenergy 
is given in terms of the in-medium and on-shell vrvr forward scattering amplitude as [B| 



oo 



S.(efc, k) = ^ J J^^if^i^p) - + ep)] / dE,^,E,^,MnAEcn.s) ■ (50) 

? n 



Here we have neglected the numerically very small contributions [^] from thermal exci- 
tations (as is usually done in the literature). The in-medium scattering amplitude, to be 
calculated from eq. (15), depends again on the pion selfenergy through the in-medium 
two-pion propagator of the intermediate state: 

G2AE, q) = L"^2Tw ■ ^51) 

eg E2-4(u;2 + E^(eg,g)) 
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Thus eqs. (8), (14), (50), (51) define a selfconsistency problem of Brueckner type wliicli 



Tlie converged results for the pion selfenergy and the in-medium M^Tr^amplitude are then 
used to calculate the thermodynamic potential and state variables as described in sect. 2 
and in sect. 3.1, respectively. According to scheme 'a' or 'b' (see sect. 2.3) we, respectively, 
either neglect or include JmS^ when evaluating the two-pion propagator eq. (51). 

3.3 Numerical Results and Discussion 

For given value of the temperature we selfconsistently calculate the in-medium vrvr am- 
plitude as well as the pion selfenergy, which are then further used to evaluate the various 
contributions to the thermodynamic potential, namely 

(i) the lowest-order skeleton graphs, eq. (12); 

(ii) the resonance-TTTT bubble skeleton graphs to all orders, eq. (28), for the resonances 
a = p(770), /2(1270) and the simulated /o(980); 

(iii) the second- and third-order skeleton graphs for p t-channel exchange, eq. (19); 

(iv) the single- (quasi-) pion contributions given by eq. (39) or (40). 

Let us first concentrate on scheme 'a' in which we neglect any imaginary part of the pion 
selfenergy. The full results for pressure-, entropy- and energy-density (full lines) are 
contrasted in fig. 1 with the results for free pions (dashed-dotted lines) and free p mesons 
(dotted lines). First we note that the quasiparticle contribution eq. (39) coincides at all 
temperatures within less than 1% with the EOS of free pions. 



The reason for that is an almost exact cancellation of the mean field contribution (second 
term in eq. (39)) and the kinematic modification caused by the replacement uJk — > 
in the exponential of the logarithm in the first term of eq. (39). But even the absolute 



we solve (at fixed temperature) by numerical iteration as discussed in ref. 0. 




(52) 
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magnitudes of these medium-induced single-particle modifications do not exceed 3% at 
any temperature. These findings are at variance with the results of ref. [0, where a 
general decrease of the thermodynamic state variables compared to the free pion gas was 
found. This decrease is due to an increase of the in-medium pion mass dominantly leading 
to a suppression of the In-term in eq. (39). The increase of the pion mass stems from 
the net repulsion of the low-energy s-wave tttt interaction (Weinberg Lagrangian). In 
our analysis, however, it turns out that it is important to include both a higher energy 
range of the vrvr interaction (essentially up to 1 GeV two-pion CMS energy) and higher 
partial waves. E.g. the resonant p-wave channel significantly contributes to -ReS^(efc, fc), 
whereas there are large cancellations in both s- and d-wave between the attractive / = 0- 



and the repulsive / = 2-channels 25 



Let us now come to the discussion of the skeleton diagrams. The effect of the first order 
graphs, eq. (12), is small: their combined contribution does not exceed 1% of Vt^J'^^iT) in 
the considered temperature range. The most important role is played by the s-channel 
p-pole graphs, see e.g. table 1. It is remarkable that nearly half of the contribution to 
PirnpiT), ST,np(T) and e.„T,p(T) is due to negative energy contributions represented by the 
C^-term in eq. (28). Without it the values for PT„rpiT), s-^T^piT) and e^7rp(^) he slightly 
below the free p-gas values p^J'^'^iT), s^p^'^'^iT) and e^p^'^^iT) (dashed-double-dotted lines in 
fig. 1). This is in agreement with the findings of ref. [|lT], i.e. that the vrvr p-wave interac- 
tion effect resembles the presence of free p-mesons, the small suppression compared to the 
free p-gas values in our case being caused by the finite width of the p-resonance (which is 
even further broadened in medium) and the missing t-channel exchange diagrams. The 
latter, when implemented in the polegraph - tttt bubble sum, would lead to a slightly 
stronger renormalization of the bare p. 

The effect of other vrvr resonances is small: pT^-^f^iT) = 3.8%pl'^'^^, pT^-^f^^T) = 0.3%pl'^'^^ 
at T=200 MeV, and similarly for entropy- and energy-density. A somewhat larger con- 
tribution is generated by skeleton diagrams from second order p t-channel exchange, see 
e.g. table 2. The third order values are already one order of magnitude smaller. 
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The overall picture emerging for the temperature dependence of the thermodynamic state 
variables is quite similar for all three: we find an appreciable enhancement of the selfcon- 
sistently calculated pressure, entropy and energy for the interacting gas compared to free 
pions for temperatures T > lOOMeV, which even exceeds the values for a free gas of pions 
and p-mesons. This excess can be traced back to negative energy contributions, naturally 
arising within our formalism, and low-energy t-channel exchange of virtual p-mesons 
between two pions. At very high temperatures (T > 200MeV) the thermal broadening 
of the TTTT resonances, in particular the dominantly contributing p(770), becomes so large 
that they no longer resemble the effect of an independent particle species (remember from 
the Beth-Uhlenbeck formalism that a sufficiently sharp two-particle resonance thermo- 
dynamically acts like an independent species corresponding to the quantum numbers of 
the resonance). 

When including a finite width, 7mE^, for the quasi-pions (= scheme 'b', as described in 
sect. 3.2.) we observe a slight overall enhancement of j9^(T), s^(T) and e^^lT) (dashed 
lines in fig. 1) compared to the results when neglecting ImT^T^. The main source for this 
increase are the second order p t-channel skeleton graphs, the p pole skeleton graphs and 
the quasiparticle contributions eq. (40), compare table 3. In case of the skeleton diagrams 
this increase is simply due to the fact that JmM,r7r (-£/') acquires a nonzero imaginary part 
below the (in-medium) two-pion threshold all the way down to £" = 0. This leads to 
an additional contribution to the £"-integral in eqs. (19), (28), accounting for the almost 
entire difference in comparison to scheme 'a'. 

We end this chapter by testing our numerical results with respect to thermodynamic 
consistency. For that we recalculate the total pressure-density by integrating the entropy- 
density according to eq. (49), which amounts to an implicit check of the stationarity 
condition, eq. (44). As can be seen from fig. 2 there is excellent agreement of the such 
calculated Pi^iT) with the 'direct' evaluation when neglecting the pion width (scheme 'a'). 
When including a finite /mE^ (scheme 'b') the values for pressure-density obtained from 
eq. (49) are slightly suppressed as compared to the 'direct' calculation, causing a loss of 
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about half the enhancement over scheme 'a'. 



4 Dropping Rho Meson Masses 

In the calculations of the previous chapter no medium effects were considered in the 
interaction kernel of the vrvr scattering equation. However, since our vrvr model is based 
on explicit (vector) meson exchange, one may expect an impact on the pseudopotentials 
due to modifications in the propagator of the exchanged mesons when exposed to finite 
temperature. The most prominent example for such a modification is the dropping of 
the vector meson masses or 'Brown-Rho Scaling' |T|, |^. Adami and Brown suggested the 
temperature dependence of the p meson mass to be 

m,(T)^m,(0) (1_^S>I\"\ (53) 



<qq>o 

where < qq >t denotes the quark condensate at finite temperature. The T dependence 
was taken as 



< qq >T=< qq >o VI - (T/Tc^)^ . (54) 

In ref. such a decrease of the p mass was implemented in a selfconsistent pion gas 
calculation within the Jiilich model as described in sect. 3.2.. Already well below 
a large accumulation of strength was found in the scalar-isoscalar channel of the tttt 
scattering amplitude close to the two-pion threshold. It is due to the strong enhancement 
of the attractive t-channel p exchange in this channel. However, it has been shown 
recently p6[ that the implementation of chiral constraints into the vrvr interaction is crucial 



to reliably calculate in-medium tttt correlations in the vicinity of the threshold. On the 
other hand, sufficiently above the two-pion threshold {Ecms > AOOMeV or so) chiral 
constraints rapidly cease to have significant impact on in-medium tttt amplitudes. This 
justifies the use of the non-chirally symmetric Jiilich model in the previous chapter, where 
no dropping vector meson masses are taken into account and thus only minor threshold 
effects in the vrvr amplitude are observed H. 
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In this chapter we will employ a chirally improved version of the Jiilich vrvr interaction |31 



2^ . It is supplemented with vrvr contact interactions as required from the gauged nonlinear 
a model |3^. To ensure the correct chiral limit for the s-wave scattering lengths a 
modified off-shell prescription for the pseudopotentials has been chosen when iterating 
the scattering eq. (14). 

Due to the additional interaction terms in the Lagrangian of the chirally improved Jiilich 
model there arise further contributions to the thermodynamic potential i7^(T). They are 
straightforwardly incorporated in the quasiparticle and in the lowest order contributions, 
fi^(T) and Q!)^1{T), respectively. For OS^^iT) we take into account all contributions up 
to second order, namely from 

- second-order p exchange as given by eq. (19) for n=2 with 

^/mM(2) = \vpex ImG^^ V,e. , (55) 

o o 

where we explicitly indicated the degeneracy factor l/4n; 

- the second-order contact interactions also given by eq. (19) for n=2, but with the 
degeneracy factor | replaced by ^ due to the absence of exchange diagrams (compare 
the remarks following eq. (4)), i.e. 

hmMi^J = ^Kont ImG^^ Ko„t ; (56) 

- the first-order p exchange plus first-order contact interaction; here the appropriate 
degeneracy factor is |: 

hmMi^J = i^ont IniG^^ V,,, . (57) 
6 6 

Numerically it turns out that higher-order contributions in p exchange or contact inter- 
actions are again negligible. 

When recalculating the EOS without any dropping of the p mass the results employing 
the chirally improved Jiilich model differ by around 5% from those presented in sect. 3.3, 
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confirming the notion that chiral symmetry has httle impact on the tttt interaction as long 
as strong threshold effects are absent. 

Let us now turn to the scenario with 'Brown-Rho Scaling' included. For simplicity we 
will neglect the coupling to the KK channel, which was shown to have a very small effect 
in the previous chapter. Relevant for our model are the decrease of the (physical) p mass 
(eq. (53)) and the corresponding reduction of the pion decay constant, 

MT) ^ rrv(T) 

A(0) m,(0) ■ ^^^^ 
The latter enters the tttt contact interactions, whereas the physical mass nip is used in 
the t-channel p exchange propagator, eq. (33). For given temperature the bare p mass 
m^p^ used in the s-channel p pole propagator, eq. (23), is adjusted such that the fully 
renormalized M^^-amplitude in free space acquires it's maximum (i.e. p resonance peak) 
at the corresponding value of mp{T) from eqs. (53), (54). For fixed mp{T), fniT) and 
W (7^) then perform selfconsistent Brueckner calculations for E^r and M^tt at several 
T-values along exactly the same fines as described in sect. 3.2. 

Tfie resulting EOS for tfie fiot interacting pion gas is displayed in fig. 3, where we have 
chosen the critical temperature of chiral symmetry restoration to be T^=YJQ MeV. For 
T>100 MeV we observe a strong increase of the thermodynamic state variables as com- 
pared to the calculations without inclusion of the Brown-Rho Scafing. Most of the effect 
is of simple kinematic origin, as can be seen by comparison to the free vr-p gas including 
the decrease of rUp. However, for temperatures below T fa 150 MeV the tttt interactions 
still generate an enhancement over the free vr-p case of up to 15%, which is quite similar to 
what we found in the pure Brueckner calculations. Just below the free 7r-p gas values 
for PniT), STr{T) and eT^{T) are close to the ones of the interacting pion gas. This is due to 
the fact that the p polegraph starts acquiring appreciable renormalization contributions 
especially from the strongly attractive t-channel exchange of a p meson of mass mp{T). 
However, in our approximation scheme for (^^^^(T), this kind of contributions is not ac- 
counted for (compare sect. 2.2.1), even though they are certainly not negligible anymore 
in the vicinity of T^. 
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5 Summary and Conclusions 

Based on the finite-temperature Greens-function formalism we have presented an analysis 
of the equation of state of a hot interacting gas of pions at zero chemical potential. 
Starting from a realistic vrvr meson-exchange model capable of describing the vacuum 
scattering data over a broad range of energies we have calculated the in-medium vrvr 
scattering amplitude and the single-pion selfenergy. The resulting selfconsistency problem 
has been solved by numerical iteration within the quasiparticle approximation (Brueckner 
scheme) 0. 

We then proceeded to evaluate the thermodynamic potential fi^(T): 

- the single-pion selfenergy was used for calculating the quasiparticle contribution 
Q^{T); the latter turned out to be very close to the values for free pions due to a 
cancellation between the mean field term and the free quasi-pion term (only when 
including a finite width for the pions a ~5% enhancement was found); 

- for the interaction part of the thermodynamic potential, Qt^t^{T), we have taken 
into account the /o(980)-, p(770)- and /2(1270)-s-channel pole graphs to all orders 
as well as the p t-channel exchange up to third order; higher orders as well as 
interference terms between s- and t-channel graphs have been estimated to be 
small and were neglected; the main contribution stems from the p pole graphs, 
their magnitude being comparable to what one expects from an admixture of free p 
mesons. 

Special attention has been paid to constraints from thermodynamic consistency, and it 
was shown that our results satisfactorily fulfill these constraints. 

In total we have found a 10-15% enhancement of ^ln{T) compared to a free gas of it and 
p mesons in the temperature range 100 MeV< T < 150 MeV. Pressure-, entropy- and 
energy-density, which were extracted from Qn{T) by means of standard thermodynamic 
relations, show a very similar behavior. This is in qualitative agreement with the re- 



sults of Welke et al. [|TT|, who employed empirical vacuum s- and p-wave tttt phase shifts 
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within a relativistic virial expansion. Our results are, however, at quahtative variance 



with the findings of Bunatian and Kampfer [|T^, who also used the finite temperature 
Greens-function approach but employed the Weinberg Lagrangian. The latter is known 
to account only for the low-energy s-wave tttt interaction, resulting in a net repulsion 
when performing the isospin weighted sum for the near-threshold tttt amplitude. As a 
consequence, their thermodynamic state variables show an overall decrease with temper- 
ature as compared to the free pion gas. 

Our conclusion from this is that the higher energy range (up to CMS energies close to 
1 GeV) as well as higher partial waves (especially the resonant p-wave) in the tttt interac- 
tion are important for a reliable description of thermodynamic properties of an interacting 
gas of pions. 

We furthermore investigated the impact of a dropping p mass according to 'Brown-Rho 
Scaling'. Using a chirally improved Jiilich model including contact interactions, arising 
from the gauged Weinberg Lagrangian, we find a similar behavior as before: even though 
the thermodynamic state variables of the interacting pion gas now exhibit a considerable 
increase near T^^, the enhancement over the free vr-p gas values (dropping p mass included) 
is not more than 10-20%. Very close to the critical temperature, where the p mass rapidly 
drops to zero, the interacting pion gas values are close to the ones of the free vr-p gas. 
This feature, however, is due to the breakdown of our approximations; in particular, 
interference terms between p s- and t-channel graphs become far from negligible, 
the strong gauge 'data', ('external') however, 
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Tables 



T [MeV] 




„ /„/ree 


■^TrTrp/ "^TT 


150 


26% 


41% 


45% 


200 


43% 


56% 


60% 



Table 1: Contributions from the p pole graph skeleton diagrams to the thermodynamic 
state variables in an interacting hot pion gas within the Brueckner calculations neglecting 
any width of the pions. 



T [MeV] 


„(2) l^free (2) 1 free (2) 1 free 
J^TTirpex/ -r TT irirpex/ it irirpex/ it 


150 
200 


7.0% 8.5% 8.9% 
9.8% 11.2% 11.7% 



Table 2: Contributions from the 2. order p t-channel exchange skeleton diagrams to the 
thermodynamic state variables of an interacting hot pion gas within the Brueckner cal- 
culations, neglecting any width of the pions. 





p%ex/pi'''' p..p/pi''' p^/pI"^'' pi^'/pfj^' 


ImE^ = 

ImT^T, < 


9.8% 43% 99.8% 158.4% 
13.3% 46.6% 107.6% 173.7% 



Table 3: Various contributions to the pressure density of an interacting hot pion gas 
at T = 200 MeV within the Brueckner calculations when neglecting (upper line) and 
including (lower line) a finite pion width. 
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Figure Captions 



Fig. 1: Thermodynamic state variables (scaled by temperature to dimensionless units) of 
an interacting hot pion gas within the Brueckner scheme with no medium modifi- 
cations applied to the two-body interaction potentials; 
upper panel: pressure-density; 
middle panel: entropy-density; 
lower panel: energy-density; 

full lines: interacting pion gas when neglecting the pion width; dashed lines: in- 
teracting pion gas with inclusion of a finite pion width; dotted lines: free p gas; 
dashed-dotted lines: free pion gas; dashed-double-dotted lines: free 7r-p gas. 

Fig. 2: Comparison of two different ways of calculating the pressure-density of an interact- 
ing pion gas within the Brueckner scheme; the full line {Iml^T^ neglected) and the 
dashed-dotted line {ImTiT^ included) correspond to the 'direct' calculations (the sum 
of eqs. (12), (19), (28) and (39)/(40)), whereas the dashed line {ImT,^^ neglected) and 
the dotted line {Itu'Etj- included) are obtained from integrating the entropy-density 
according to eq. (49) . 

Fig. 3: Thermodynamic state variables (scaled by temperature to dimensionless units) of 
an interacting hot pion gas within the Brueckner scheme plus additional medium 
modifications of the two-body pseudopotentials in form of dropping mp{T) and 
/7r(T') ; here, a chirally improved version of the Jiilich tttt interaction has been 
employed; 

upper panel: pressure-density; 
middle panel: entropy-density; 
lower panel: energy-density; 
line identification as in fig. 1. 
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